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Much new solid state technology for single-photon sources 1 , detectors 2 3 , photovoltaics 4 
and quantum computation 5 6 relies on the fabrication of strained semiconductor nanos- 
tructures. Successful development of these devices depends strongly on techniques 
allowing structural analysis on the nanometer scale. However, commonly used mi- 
croscopy methods 7 10 are destructive, leading to the loss of the important link between 
the obtained structural information and the electronic and optical properties of the 
device. Alternative non-invasive techniques such as optically detected nuclear mag- 
netic resonance (ODNMR) so far proved difficult in semiconductor nano-structures 
due to significant strain-induced quadrupole broadening of the NMR spectra. Here, 
we develop new high sensitivity techniques that move ODNMR to a new regime, 
allowing high resolution spectroscopy of as few as 10 5 quadrupole nuclear spins. By 
applying these techniques to individual strained self-assembled quantum dots, we 
measure strain distribution and chemical composition in the volume occupied by the 
confined electron. Furthermore, strain-induced spectral broadening is found to lead 
to suppression of nuclear spin magnetization fluctuations thus extending spin coher- 
ence times. The new ODNMR methods have potential to be applied for non-invasive 
investigations of a wide range of materials beyond single nano-structures, as well as 
address the task of understanding and control of nuclear spins on the nanoscale, one 
of the central problems in quantum information processing 1 ^^. 

Most nuclei used in optically active III-V semiconductor nanostructures possess non-zero 
quadrupole moments sensitive to electric field gradients caused e.g. by strain. Due to strong 
spectral broadening NMR on quadrupole nuclei is challenging even for macroscopic samples"^. 
As a result control of nuclear spins using NMR has only been achieved in strain- free GaAs/AlGaAs 
semiconductor quantum dots (QDs) — ~—. By contrast application of similar techniques to widely 
researched strained self-assembled quantum dots were limited to large ensembles of QDs, where 
quadrupole broadening lead to uncertainty in interpretation of the effects observed under radio- 
frequency (rf) excitatio n 20 ' 21 . 

The new ODNMR spectroscopy technique we present is based on continuous-wave broadband 
rf excitation with specially designed spectral patterns, which for nuclei with spin / in strained 
structures provides sensitivity enhancement by a factor of (J + 1/2) 3 compared to conventional 
saturation NMR techniques previously used in unstrained nano-structures^^. Such enhance- 
ment (particularly large for high spin nuclei such as 1=9/2 indium) brings NMR to a qualitatively 
new level: it allows high resolution spectroscopy in individual few-nanometer-sized strained quan- 
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turn dots having only 10 4 -10 5 nuclear spins which are subject to large inhomogeneous quadrupole 
broadening up to 20 MHz. To demonstrate the capabilities of this technique we use NMR spectra 
measured for InP/GalnP and InGaAs/GaAs QDs to estimate several key structural properties 
such as material composition, elastic strain magnitude and distribution in the volume sampled by 
the electron wave- function. Furthermore, using related techniques we find that spectral broaden- 
ing caused by strain results in enhancement of the coherence of nuclear spins by a factor of 10 
compared with unstrained structures, leading to coherence times of ~ 2.5 ms. 

In what follows we present results for two different types of strained semiconductor nano- 
structures: InP/GalnP and InGaAs/GaAs self-assembled quantum dots [for sample details see 
Supplementary Information (SI) Section SI]. All measurements were performed at T = 4.2 K, in 
external magnetic field B z normal to the sample surface. Under excitation with circularly polarized 
light, nuclear spins become strongly polarized due to spin transfer from electrons via the hyperfine 
interaction- 2 . The resulting nuclear spin polarization on the dot is detected in photoluminescence 
(PL) of excitons in single QDs as shown in Fig. []Ja) for InGaAs and in Fig. DJb) for InP dots 
in high magnetic field B z > 5 T. Each spectrum consists of an exciton Zeeman doublet with 
splitting E z . Detection of changes in E z allows measurement of the electron Overhauser shift^ 1 ^ 
proportional to the degree of nuclear polarization P/y. ODNMR measurements were carried out 
using the pump-probe method schematically shown in Fig. [3Jc). The dot is first excited with a 
circularly polarized laser pulse of duration T pump . After that the rf excitation is applied in the 
dark for duration T r f. Finally a short (T pro b e ) laser pulse is applied to measure the PL spectrum 
and probe the effect of the rf field on nuclear polarization, which allows NMR signal to be obtained 
as the absolute magnitude of the Overhauser shift. More details of experimental methods can be 
found in SI Sec. S2. 

The explanation of our method and its comparison with the "saturation" NMR techniques 
applied in strain- free material o 17 ' 18 is given in Fig. [21 using an example of spin 1 = 5/2 nuclei. 
In an external magnetic field B z along the Oz axis nuclear spin levels experience Zeeman shifts 
oc B Z I Z determined by their spin projections I z . The oscillating magnetic field B r f perpendicular 
to Oz couples only spin levels with I z differing by ±1. If the nucleus is subject to an electric 
field gradient (EFG) along the Oz axis (e.g. induced by elastic strain), in addition to the splitting 
induced by B Z: the energies of spin levels will change 8 - 1 ^ by a value proportional to ij (see SI Sec. 
S3 ), and all dipole active transitions will have different frequencies, as depicted in Fig. [2jb). 

Let us consider an ensemble of nuclei with spins I = 5/2 all subject to the same EFG and 
polarized, so that I z = 5/2 levels have higher population than I z = —5/2. The probabilities pi z 
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FIG. 1. (a) and (b) show typical photoluminescence spectra for InGaAs/GaAs (a) and InP/GalnP (b) 
quantum dots measured in magnetic field B z — 5.3 and 6 T, respectively. Squares (circles) show PL spectra 
measured for er~( + ) excitation exhibiting differing Zeeman splittings E z due to the nuclear spin polarization 
(anti-) parallel to B z . (c) shows the timing diagram of the experimental cycle including optical pump and 
probe, and rf excitation. We use T pump = 4 7 s, T r f = 4 6 s and T pro fc e = 3 -f- 16 ms. 



to find a nucleus with spin I z will depend on I z as sketched by the solid lines in Fig. [21(c). The 
total nuclear spin polarization degree (which is detected optically) is Pjy = J2k=-iPk x k/I, so 
that \Pn\ < 1- Arrows in Fig. [2(c) indicate the maxima in the NMR spectrum in Fig. [2](b) 
corresponding to allowed transitions between adjacent pairs of nuclear spin levels. 

In conventional "saturation" NMR spectroscopy, radio- frequency excitation at a frequency v 
or, a distribution of frequencies with a width w exc is applied [Fig. [2(a)] . Pn will only change 
in the case when v is in resonance with a transition between I' z and I' z + 1 levels, for example, 
—3/2 <-)■ —1/2 in Figs. [2](a-c) . This occurs as the populations of these spin levels equalize under 
sufficiently long resonant rf excitation [dotted lines in Fig. [2(c)], which at the same time has no 
effect on populations of the other spin levels. As a result, the overall change in nuclear polarization 
is small, making the resonance difficult to detect. 

Very major enhancement of the changes in Pjy can be achieved by using an alternative approach 
developed in this work [see Fig. L2(d)]: we use broad band excitation with a continuum spectrum 
containing a gap of width w gap . As the rf excitation spectrum in Fig. [2(d) is an inversion of that 
in Fig. [2(a), we introduce the term "inverse" spectroscopy. The effect of such excitation on the 
populations pi z is demonstrated in Fig. [2(f). If the gap is out of resonance with all transitions, 
all pj z are equalized (solid lines) and nuclear spin polarization is completely erased (Pjv = 0). If, 
however, the gap is in resonance with the V z I' z + 1 transition, i. e. one of the transitions 
is not excited, the equalization of populations occurs independently for two groups of levels with 
I z < I' z or I z > I' z + 1 (dotted lines). Thus the resonance condition for one of the transitions affects 
the populations of all states including I z = ±1 states, which give the largest contribution to Pn- 
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FIG. 2. Schematics explaining "saturation" (a)-(c) and "inverse" (d)-(f) NMR spectroscopy applied to 
quadrupolar nuclei. The case of an ensemble of spin-5/2 nuclei subject to the same electric field gradient 
is considered for clarity, (a) and (d) show rf excitation spectra in the two spectroscopy methods. Spectra 
of dipole transitions between the nuclear spin levels are shown in (b) and (e). (c) and (f) show population 
probabilities of the nuclear spin levels. Arrows show the transitions in the nuclear spectra in (b) and (e) 
corresponding to pairs of the spin states coupled by the rf field. Solid lines show the population probabilities 
for the case when the rf maximum in (a) and the gap in (d) are off resonance with all transitions. Dashed 
lines show the case when the maximum in (a) and the gap in (d) are in resonance with —3/2 ■<-> —1/2 
transition. 



In the experiment, the "inverse" NMR spectrum is obtained by scanning the central frequency 
v of the gap, while w gap is chosen to control the balance between the spectral resolution and 
sensitivity (NMR signal amplitude). It is possible to show (see SI Sec. S4 ) that for "inverse" 
method the enhancement of the changes in nuclear polarization exceed (I + 1/2) 3 (=125 for spin 
/ = 9/2) compared to the saturation NMR method in Fig. E)^a). This is a significant improvement 
compared to the existing "population transfer" technique where a maximum enhancement of 21 
can be achieved 16 . 

Fig. [3] shows a set of ODNMR spectra measured on single InP (a,b) and InGaAs (c-f ) QDs at 
B z ~ 5.3 T. In the nominally InP dots grown in GalnP barriers, the "inverse" technique allows 
to resolve contributions from quadrupole nuclei 115 In, 69 Ga and 71 Ga within the volume probed 
by a single electron. In Fig. [3{a) the 115 In peak dominating the spectrum, consists of a sharp 
central line [corresponding to —1/2 f-> 1/2 central transition (CT)] with amplitude ~40 ^eV at 
~49.7 MHz, and two broad bands of satellite transitions (STs) to lower and higher frequencies each 
stretching up to ±7 MHz. These sidebands are due to strain-sensitive, quadrupole split transitions 
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FIG. 3. ODNMR spectra measured at B z ss5.3 T in self-assembled QDs using the "inverse" method with 
<j + (red curves) or er~ (blue curves) optical pumping and using the "saturation" method with rj + pumping 
(black curves), (a) and (b) show ODNMR spectra for InP/GalnP dots measured with resolution w gap =280 
kHz in (a) and 16 kHz in (b), and applying "saturation" spectroscopy with a width w exc =450 kHz for 115 In 
and single-frequency excitation for spin-1/2 31 P in (a), (c)-(f) show ODNMR spectra for InGaAs/GaAs 
dots measured with resolution w gap =280 kHz in (c), 8-24 kHz in (d), and 800 kHz in (e,f). Vertical line in 
(b) shows indium frequency vi n 49633 kHz corresponding to unstrained InP. 



between spin levels with \I Z \ > 1/2. The relative amplitudes of the side-bands reflect the initial 
alignment of the nuclear spins by circularly polarized excitation: the high (low) frequency band 
has a higher (lower) intensity for o~ (cr + ) excitation. By contrast, with the saturation method 
applied to 115 In nuclei a weak line with amplitude ~3 //eV [see Fig. [3(a)] can be detected only 
for rf excitation width w exc = 450 kHz or larger, i.e. at low resolution. For spin-1/2 31 P nuclei 
unaffected by quadrupole effects the saturation method with monochromatic rf excitation reveals 
a single line with a width of ~8 kHz at vp ~91605 kHz. 

The ODNMR spectra of a InGaAs/GaAs QD measured using the "inverse" method (u> ffap =280 
kHz) are shown in Fig. [3(c). Here the central transitions have similar amplitudes for the four 
isotopes present in the dot, revealing significant substitution of indium by gallium. Satellite transi- 
tions are observed only for spin-9/2 115 In and have similar widths to the case of InP dots. In order 
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to observe these transitions for spin-3/2 nuclei we have carried out "inverse" NMR measurements 
with a larger w gap =800 kHz as shown in Figs. |3](e) and (f) for 71 Ga and 75 As respectively. For 
71 Ga the spectral range where STs are observed is within ±2.5 MHz on both sides of the central 
line, whereas it is significantly broader for 75 As due to its larger quadrupole moment Q. 

The structure of the central transitions is affected by strain only via weak second-order 
quadrupole interaction-. It is measured using "inverse" NMR with a smaller w 9ap =8H-24 kHz 
to provide higher resolution. The smallest linewidth of ~8 kHz is observed for 71 Ga in InGaAs 
QDs while more sensitive to quadrupole interaction 75 As as well as 115 In in both materials have 
linewidths of ~40 kHz. Broadening of the CT line as well as its shift with respect to resonance 
frequency in unstrained bulk material [shown by vertical line in Fig. E^b) for 115 In in InP] appear if 
electric field gradient direction deviates from that of magnetic field and thus reflect the distribution 
of elastic strain directions within the volume probed by the electron. 

Quadrupole effects are insensitive to isotropic (hydrostatic) strain, but can provide information 
on biaxial and shear strains. For a given electron wavefunction and strain distribution within the 
QD it is possible to calculate an "inverse" NMR spectrum that can be directly compared with 
experiment. Since detailed quantum dot modeling^ is beyond the scope of this work we limit our 
discussion to some quantitative estimates. For example, it is observed in Figs. [3](a,c) that satellite 
transitions I' z I' z + 1 of 115 In corresponding to different V z are not resolved. This signifies a strong 
variation of the quadrupole shifts over the volume of the dot resulting from the variation of elastic 
strain. In particular, it can be seen that ST bands have non-zero amplitudes at the CT frequency, 
implying that for some nuclei the quadrupole splitting is zero. This can be explained by relaxation 
of the biaxial strain near the center of the dot as predicted by model calculations^ 1 ^. On the 
other hand, the maximum values of the biaxial strain can be readily estimated from the maximum 
frequency shifts observed for the STs. For the spin-3/2 71 Ga, the largest shift of the STs from the 
CT in InGaAs/GaAs QDs is ~2.5 MHz [Fig. ffle)]. Using the values of the EFG-elastic tensor 
measured in bulk GaAs^ this allows the maximum biaxial strain with principle axis along Oz to 
be estimated as |e&| ~ 6%. This is somewhat smaller than the maximum \eb\ ~ 9 -f- 15% predicted 
for interfacial regions of InAs dots of different shapes^, possibly a signature of "smoothing" of the 
interface between the dot and the barrier due to interdiffusion. Using NMR on 115 In we find in a 
similar way the maximum strain for InP/GalnP dots of | e^, | ~ 5%. (See details in SI Sec. S6 ). 

The NMR data in Fig. [3] can be used to estimate relative gallium and indium concentrations 
PGa an d pin- For InGaAs dots in Fig. EJc) similar CT signals are observed for Ga and 115 In 
isotopes. However, we need to take into account, that for the 115 In isotope the NMR signal is 
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enhanced due to the larger nuclear magnetic moment of indium /inis j n \ / fJ>aiG a ) ~ 2, and its larger 
natural abundance (~96%) compared to 40% for 71 Ga. In this way we estimate pj n 20% and 
PGa ~ 80% for InGaAs dots. Similar estimates for InP dots give pj n « 65% and pc a ~ 35%. These 
estimates are further confirmed by more detailed measurements and analysis presented in SI Sec. 
S5 . 

The spectral widths exceeding 10 MHz observed for spectra in Fig. [3] are due to inhomogeneous 
quadrupole broadening. However, each nuclear spin transition has a finite intrinsic linewidth w nuc 
determined by fluctuations of the local fields, created by neighboring nuclear spins. In unstrained 
structures w nuc can be determined by spin-echo techniques^, using ir and ir/2 rf-pulses for coherent 
manipulation of nuclear spins. In strained structures uniform excitation of nuclear spins, with 
spectral dispersion of several MHz, would require very short rf-pulses, which in turn require rf- 
powers that are difficult to achieve with available hardware^. However, as we now show, nuclear 
spin coherence can be probed using an alternative approach based on low-power continuous-wave 
rf excitation with a special spectral pattern. 

This approach is demonstrated in Fig. H]for InP QDs. We measure the rate R r f of nuclear spin 
depolarization induced by broad-band rf excitation consisting of equally spaced spectral modes 
[rf spectra are sketched with solid lines in Figs. |4] (b,c)]. The measurements are carried out for 
different spacings w m between the modes, while keeping a constant spectral range 44-47.5 MHz 
corresponding to the satellite transitions of In nuclei [see Fig. Eta)] and constant total power of 
the rf excitation. 

The strong reduction of R r f with increasing w m observed in the experiment [Fig. Ufa)] is 
qualitatively explained in Figs. EIb,c). For w m < w nuc [Fig. [31(b)] the discrete structure of the 
rf-spectrum is averaged out by the broadening of the nuclear spin transition and each nuclear 
spin transition is excited by many modes resulting in fast depolarization. If the mode spacing 
is increased (w m > w nuc ) there is a high probability that a given nuclear spin transition is out 
of resonance with all modes [Fig. Etc)]. Strong suppression of R r f observed experimentally at 
w m ~4 kHz [Fig. [U^a)] will take place for very large mode spacing w m ^> w nuc when most of the 
nuclear transitions will be out of resonance with rf excitation. Model fitting allows an intrinsic 
linewidth w nuc ~0.13 kHz to be determined. Furthermore, analysis of the lineshape of individual 
nuclear spin transitions reveals a deviation from a Lorenzian lineshape indicative of non-exponential 
dynamics of the transverse nuclear magnetization (see details in SI Sec. S7 ). 

From w nuc we can estimate the nuclear spin decoherence time of indium nuclei as Ti ~ 
l/(iTw nuc ) 2.5 ms. This exceeds by nearly an order of magnitude previously reported T2 for 
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FIG. 4. (a) Experimental dependence of rf induced nuclear depolarization rate R r f on the spacing w m 
between spectral modes of the rf excitation, (b) and (c) qualitatively explain the effect of the variation of 
w m on the nuclear polarization decay rate (see text). Solid and dashed lines show schematically spectra of 
rf excitation and nuclear spin dipole transition, respectively. 



quadrupole nuclei in unstrained GaAs/AlGaAs QWs^S (0.27 ms) and QDsi^ (0.35 ms). This can 
be interpreted as a result of large quadrupole shifts, caused by inhomogeneous elastic strain, that 
make flips of spin states with \I Z \ > 1/2 energetically forbidden, partially freezing fluctuations of 
the dipole-dipole field. The nuclear spin decoherence time is likely to be determined by spin flips 
of the —1/2 -H- 1/2 transition of both quadrupole nuclei and phosphorus. Thus, while leading to 
inhomogeneous spectral broadening exceeding 10 MHz, quadrupole interaction has the beneficial 
effect of narrowing of individual transitions down to w nuc «0.13 kHz. 

The proposed techniques may become a useful tool in development of single QDs for exper- 
iments on electron and hole spin coherence, as direct measurements of material properties such 
as composition and strain have now been made possible. This will potentially allow a generation 
of semiconductor nano-structures with tailored material properties for these applications. Fur- 
thermore, our experiments reveal very robust coherence properties of nuclear spins in strained 
structures, showing that introduction of strain is one of the ways to achieve a very stable nuclear 
spin bath, which is of potential significance for enhancing electron and hole coherence on the 
nanoscale. Finally, we note that these techniques may be applicable in a wide range of NMR 
detection schemes and structures beyond quantum dots. One of the challenges will be to develop 
nano-NMR techniques (including imaging) sensitive to yet smaller numbers of nuclei in the solid 
state environment. 
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SUPPLEMENTARY INFORMATION 

In this document we refer to the figures 1-4 of the main text as well as to supplementary figures 
S1-S5. 

The document consists of the following sections: 
IS 11 Details of sample structure and growth, 
IS21 Details of experimental techniques, 
IS3L Nuclear spin spectrum in the presence of strain, 
IS41 Calculation of the NMR signal for "inverse" spectroscopy, 
IS51 Estimation of the chemical composition of the dots, 
IS61 Estimation of strain in a QD, 

IS71 Spin coherence in the nuclear spin ensemble subject to inhomogeneous strain. 

SI. DETAILS OF SAMPLE STRUCTURE AND GROWTH 

We use InP/GalnP self-assembled quantum dots (QDs) grown by metal-organic vapor-phase 
epitaxy (MOVPE) and InGaAs/GaAs quantum dots grown by molecular beam epitaxy (MBE). 
Both samples are not intentionally doped and have no electric gates. 

The InP/GalnP sampl o ' was grown in a horizontal flow quartz reactor using low-pressure 
MOVPE on (100) GaAs substrates misoriented by 3° towards (111). The growth temperature of 
the GaAs buffer and bottom Gao.5Ino.5P layer was 700° C. Before proceeding to the deposition of 
InP and the Gao.5Ino.5P capping layer, the wafer was cooled to 650° C. The grown GalnP layers 
were nominally lattice matched to GaAs. A low InP growth rate of l.lA/s and deposition time of 
10 seconds were chosen. 

The InGaAs/GaAs sample^— consists of a single layer of InAs quantum dots (QDs) placed 
within a microcavity structure which is used to select and enhance the photoluminescence from part 
of the inhomogeneous distribution of QD energies. The sample has been produced using a single 
step MBE process and consists of a GaAs cavity of thickness A/n formed between an asymmetric 
set of distributed Bragg reflector pairs, which uses 16 pairs of GaAs/Alo.sGao^As distributed Bragg 
reflector pairs below and 6 pairs above the cavity. The resulting cavity Q factor is ~250 and the 
cavity has a low temperature resonant wavelength at around 920 nm. The luminescence of the 
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QDs is further enhanced by the presence of a GaAs/AlAs short-period superlattice surrounding 
the QD layer. The superlattice and DBR layers are grown around 620° C. The quantum dots were 
formed by deposition of 1.85 monolayers (MLs) of InAs at a growth temperature of 510° C and a 
growth rate of 0.1 ML/s. The deposition amount is just above that required for the nucleation of 
QDs (~1.65 ML) but is well below the amount required to produce a mature QD distribution. As 
a result, we obtain a low density of infant QDs at the post-nucleation stage, which are small and 
have a low concentration of indium. 

Both QD samples were examined using transmission electron microscopy (TEM) . Fig. [Si] shows 
images for InP/GalnP (a-b) and InGaAs/GaAs (c-e) samples taken under dark field 002 condition, 
which gives compositionally-sensitive diffraction contras1>£&. We find that InP dots are disk shaped 
with lateral size of ~75 nm and a height of ~4.5 nm. Images taken on the InGaAs/GaAs sample 
clearly show the AlGaAs layers of the Bragg mirror and the superlattice as well as pyramidal 
shaped QDs approximately 30 nm wide at their base and 5 nm high in the center. 



S2. DETAILS OF EXPERIMENTAL TECHNIQUES 



The experiments are performed with the sample placed in an exchange-gas cryostat at T = 4.2 K, 
and using an external magnetic field B z normal to the sample surface. In order to detect nuclear 
polarization on the dot we use high resolution micro-photoluminescence (PL) spectroscopy of single 
QDs (see experimental setup scheme in Fig. IS2|h The QD PL is excited by a laser resonant with 
the wetting layer states (E exc =1.88 eV for InP dots and E exc =lA6 eV for InGaAs dots) and 
analyzed with a 1 meter double spectrometer coupled to a CCD. Manipulation and probing of 




FIG. SI. Transmission electron microscope (TEM) images showing the structure of InP/GalnP (a-b) and 
InGaAs/GaAs (c-e) QDs. (b) shows zoomed part of the image in (a), (c) and (e) are zoomed parts of (d). 



3 




Arbitrary 
waveform 
generator 



rf power 
amplifier 



Liquid He 2 bath 



FIG. S2. Experimental setup. The sample is placed inside an insert filled with low-pressure helium gas with 
optical access from the top. The insert is immersed in a liquid helium bath cryostat. Photoluminescence 
(PL) is excited by diode lasers. All experiments are performed in the Faraday geometry with excitation and 
collection normal to the sample surface. Focusing of the laser and collection of PL from a QD are performed 
with an aspheric lens with focal length / ~ 2 mm and numerical aperture NA« 0.5. The collected PL signal 
is coupled via an optical fiber to the entrance slit of a double spectrometer with a charged coupled device 
(CCD) camera. Using line fitting or calculating center of mass it is possible to determine the shifts in QD PL 
energy as small as 1 /ieV. We use two identical lasers combined via a beamsplitter which allow independent 
control of power and polarization of "pump" and "probe" pulses. Mechanical shutters are used to control 
the light of both lasers and PL signal according to the pump-probe timing diagram [Fig. [IJc) of the main 
text]. Time accuracy of the shutters is ~1 ms. Magnetic field up to 8 T can be applied perpendicular to 
the sample surface (along the Oz direction). Oscillating radio- frequency (rf) magnetic field along the Ox 
direction is produced by two coils on each side of the sample. The rf signal is produced by digital arbitrary 
waveform generators with analog modulators, and after the power amplifier, coupled to the coil via a coaxial 
cable. 

nuclear spin polarization relies on the hyperfine interaction of electrons and nuclear spins"' 7 , and 
requires polarization-resolved excitation and detection of light as described in the main text. 

The waveform for the radio-frequency excitation is produced by a digital arbitrary waveform 
generator. "Inverse" spectroscopy with a continuous rf excitation spectrum shown in Fig. [21(d) 
requires an aperiodic signal, which can not be produced by a digital device. We thus approximate 
the spectrum shown in Fig. L2Jd) with a spectrum consisting of a large number of discrete modes 
with equal spacing w m (typically w m = 0.4 kHz) as shown in Fig. IS3I Since the rf signal has a finite 
power, it has to be limited in the spectral domain to a width w exc (typically w exc = 20 MHz). The 
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FIG. S3. Schematics of the spectrum of rf excitation. The spectrum has a total width w exc up to 20 
MHz and consists of many modes with equal amplitudes B m and spacing w m , which is varied in different 
measurements in the range 0.1 -j- 4 kHz (w m = 0.4 kHz is used for "inverse" spectroscopy). For "inverse" 
NMR spectroscopy the spectrum can also have a gap in the center with a width w gap varied in the range 
8-800 kHz. 



signal waveform is synthesized to have no spectral components within a gap of width w gap . However 
due to imperfections of the rf circuit (harmonics and spurious noise) the amplitude of the rf field 
within the gap is not strictly zero. In our experiments the spectral density of the rf power inside 
the gap is ~1000 times smaller than the spectral power density of the modes. The typical mean 
square amplitude of the in-plane rf oscillating magnetic field used in the "inverse" spectroscopy 
experiments is yj {B^) ~ 0.15 mT, while the phases of the modes are chosen to minimize the crest 
factor so that the peak value is B r j < 0.3 mT. In our experimental setup continuous rf excitation 
of such amplitude results in sample heating of less than 1 K. 

The duration of the rf pulse for "inverse" spectroscopy T r f = 5.5 s is chosen to be long enough 
to produce the steady-state population probability distribution of the nuclear spin states shown in 
Fig. [2](f) ■ On the other hand for spectroscopy on phosphorus we use T r y=50 ms which is shorter 
than the time required to completely erase nuclear polarization, and thus gives an unsaturated 
spectrum allowing broadening to be avoided. Both InP and InGaAs dots used in this work exhibit 
long nuclear spin decay tiniest T\ >100 s so that natural decay of nuclear polarization during rf 
excitation and probe pulse is negligible. 

The "inverse" spectra (see Fig. [3]) measured with a + optical pumping were calculated as 
a difference E 9 z av — E™°9 a P Q f the exciton Zeeman splitting E 9 z ap measured using rf excitation 
with a gap in the spectrum at frequency v and the splitting E™° ' 9ap measured without the gap. 
Such correction allows NMR signal to be expressed as an absolute value of the hyperfine spectral 
shift. For a~ pumping the difference was taken with the opposite sign (£f ' gap — E§ ap ) to simplify 
comparison with spectra measured under a + pumping. For "saturation" spectroscopy cr + pumping 
was used and the signal was calculated as E™° r ^ — E^ , where Ez(E%° is the splitting measured 
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FIG. S4. Diagram of nuclear spin energy levels for the case of 7=3/2 spin. Dipole transitions are shown 
with arrows. Zeeman interaction leads to equidistant shift of spin levels resulting in a degenerate transition 
at a Larmor frequency vl ■ Quadrupole interaction, which can be treated as a perturbation, leads to further 
shifts. The first order correction makes transition frequency vl + v\i dependent on the spins of the 
levels I' z and I' z + 1 thus removing the degeneracy. However the correction v_*j 2 1 / 2 to the frequency of the 
central transition (transition between -1/2 and 1/2 states) appears only in the second order. 

with (without) rf excitation centered at frequency v. 



S3. NUCLEAR SPIN SPECTRUM IN THE PRESENCE OF STRAIN 

Below we briefly summarize the effect of external fields on the nuclear spin spectrum. In 
magnetic field B z along the Oz axis and in the presence of elastic strain, the Hamiltonian for a 
nuclear spin I can be written as^: 

Hnuc = ~hv L I z + H Q , (SI) 

where v L = r yB z /(2ir) is Larmor frequency, h - Planck constant, 7 is nuclear gyromagnetic ratio 
and Hq describes interaction of the nuclear quadrupole moment with the electric field gradient 
(EFG), described by a second rank traceless tensor of the electrostatic potential second derivatives 
Vij. In the frame Ox'y'z' with the axes along the principal axes of Vif 

H Q = hu Q [3/ 2 2 , -I 2 - V (I 2 X , - I 2 )] /6, (S2) 

where vq = 2l{2l-\)h anc ^ ^ = Vx ' x y ; V f v> describe strength and deviation of the EFG from the 
axial symmetry, respectively 38 . 

In the case of high magnetic field (ui S> vq) studied here the quadrupole interaction can be 
treated as a perturbation. The effect of external fields on nuclear spin level energies for the case 
of spin 1=3/2 is shown in Fig. IS41 Without quadrupole interaction frequencies of all transitions 
equal to vl determined by Zeeman energy. In the first order of vq, the frequency of the I' z I' z + 1 
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transition, vl + Vp T ,,-,, becomes dependent on I' . For the uniaxial EFG with the main axis 
along the external magnetic field this shift reads as vpp+i = u q{I' z + 1/2) ; i-e. 2/ NMR lines are 
observed equally spaced by vq. In the case of non-axial symmetry of the EFG (77 7^ 0) or non-zero 
angle 9 between EFG principle axis Oz' and magnetic field direction Oz, further changes in the 
transition frequencies are observed and the first order shift reads as: 

vjpj, +1 = vq{I' z + 1/2) (1 +3 cos 2(9 - 2ry cos 20 sin 2 0)/4, (S3) 

where eft is the angle describing orientation of EFG axes with respect to Ozz' plane. For example, 
it follows from Eq. IS3I that for 77 = the first-order quadrupole shift can be canceled if strain axis 
deviates from the magnetic field by an angle 9 ~ 54.7°. 

However, the shift of the central transition —1/2 o 1/2 appears only in the second order of vq 
and can be written as: 

2 

^-1/2,1/2 = 9^( J ( / + 1 ) " 3/4)G(M,0), (S4) 

where G{9 = 0, 77 = 0, <p) = and — 1 < G(9, 77, 0) < 1/2 for all possible values of 9, 77, [complete 
expressions for G(6,r),<f>) can be found in Ref.— ]. Thus for ul 3> the shift of the CT is much 
smaller than for the STs, resulting in a narrow central peak in the NMR spectra. 



S4. CALCULATION OF THE NMR SIGNAL FOR "INVERSE" SPECTROSCOPY 

Let us first consider "saturation" spectroscopy on quadrupolar nuclei subject to EFG when the 
spectrum of rf excitation has a single frequency component only. If the rf frequency v is equal to 
the frequency of the transition between I' z and I' z + \ states [—3/2 <-> —1/2 in Figs. El^a-c)], the 
populations of these spin levels are equalized [the dotted lines in Fig. EJ^c)]. The populations of 
all other spin levels are not affected, and the total change of nuclear spin polarization Pn due to 
rf excitation is determined only by the contributions of I' z and I' z + 1 states. Taking the difference 
between the initial polarization [(I' z + l)pj/_|_i + I' z pp]/I and the polarization after saturation of 
the transition [(I' z + 1) + I' z ] x (pi' z +i +pj/)/(2I), we find the amplitude of the detected NMR signal 
as: 

AP^ = (p Ii+1 -p r J/(2I). (S5) 

By contrast in the "inverse" spectroscopy method [Figs. [2^d)-(f)] all spin states contribute to the 
NMR signal. For simplicity we consider the case of large EFG and narrow gap in the rf spectrum, 
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i. e. when the frequency of only one nuclear transition coincides with the gap. Calculations for a 
general case are bulky but straightforward. First we note that if the gap is not in resonance with 
any transitions then all nuclear polarization will be erased yielding the final polarization P/v = 0. 
If the gap is in resonance with the I' x •<-)■ I' z + 1 transition [—3/2 -H- —1/2 in Figs. [2] (d-f)], i. e. 
this transition is not excited, then the equalization of populations takes place separately for two 
groups of spin levels: for states with I z < I' z and for states I z > I' z + 1 [dotted lines in Fig. E^f)]. 
This is because transfer of population induced by the rf field is allowed only between spin levels 
with I z differing by ±1. For each group, the population probabilities of the spin levels after the 
rf pulse will be the average of the their initial populations pj z . Thus non-zero polarization AP™ V 
will be retained: 

Apr = ^E»+ L± ^ ±1 £ p»- (SB) 

k=-I k=I' z +l 
The initial polarization of the nuclear spins created via the hyperfine interaction with the opti- 
cally pumped spin-polarized electrons can be expressed in terms of the electron spin temperature 
T e , assuming that the nuclear spin populations pi z follow a Boltzman distribution. In large mag- 
netic fields when vl 3> uq, the energies of the nuclear spin levels are approximately proportional 
to I z and the initial populations can be described as*^: 

Ph = Z ~ X exp(7 2 /3), 

where Z is the normalization factor, ks - Boltzman constant, E e z - the electron Zeeman splitting. 
f3 describes the dynamic nuclear polarization occurring as a result of deviation of the electron spin 
temperature T e from the bath temperature T due to optical orientation of the electrons. In the 
studied case |T e | <C T and the term E e z / '(ksT) can be neglected. 

Using the Boltzman population distribution of Eq. IS7I we can calculate the sums in Eqs. IS5I 
and IS6I and obtain the following expressions for the signal amplitudes for both "saturation" and 
"inverse" NMR: 

ApSa t,B = [e (J - +1)/? ~ e 1 '^} smhQ3/2) 
N 21 sinh[(2/ + l)/3/2]' 

a p inv,B _ (I ~ 4) + {I + I' Z + l)e (2J+1)/3 ~ (21 + 

~ 2/ [ e (2/ + l )/3 _ J] • 

It is useful to consider several practical cases. For example, the amplitude of the —1/2 o 1/2 
central transition (CT) signal is obtained from Eqs. [S8]by substituting I' z = —1/2. For a given 
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/, the CT NMR signal APjv becomes a function of f3 and it can be easily derived that for any /3, 
"inverse" spectroscopy yields a signal enhancement of at least APt^ v,B / APf^*'" 8 > (J + 1/2) 3 . For 
example, for 1 = 9/2 spin (indium) it is AP^/AP^ > 125. Since the CT exhibits a narrow 
line in the NMR spectrum it is much easier to detect even in the presence of large quadrupole 
effects. Thus the use of "inverse" NMR greatly enhances the sensitivity for detection of even 
small amounts of quadrupole isotopes, in particular with large I. For the smallest half-integer 
quadrupole nuclear spin I = 3/2, "inverse" NMR for the CT still leads to a signal enhancement of 
AP™' S /AP^ af,B > 8. For the satellite transitions (STs) the sensitivity of "inverse" NMR reduces 
with increasing I' z . However, even for the most split-off transitions —I f-) —I + 1 and / -If)/, 
the enhancement is AP™ V,B / AP™ t,B > 9 for / = 9/2 and > 3 for / = 3/2. 

It is also useful to compare NMR sensitivity in strained and unstrained structures. In the 
absence of quadrupole effects all nuclear transitions have the same frequencies, and thus for both 
"saturation" and "inverse" techniques the NMR signal will be given by the initial polarization 
degree Pjv- Using Eqs. IS7I and IS8I we can calculate the ratio of the CT NMR signal Apy in 
the presence of EFG and the total signal P/v for zero EFG. For "inverse" NMR we find that for 
all this ratio is AP™ V ' B /P N > 0.55 for / = 9/2 (AP^ V ' B /P N > 0.66 for / = 3/2), implying 
that sensitivity in the strained structures is only two times smaller compared to the structures 
with zero quadrupole effect. By contrast, for saturation spectroscopy, non-zero EFG results in 
suppression of NMR signal particularly strong for nuclei with large spin: we find that at least 
AP s N at ' B /P N < 0.0061 for / = 9/2 and AP^ at ' B /P N < 0.1 for / = 3/2. 



S5. ESTIMATION OF THE CHEMICAL COMPOSITION OF THE DOTS 

This section presents an experimental method and numerical analysis enabling estimation of Ga 
and In intermixing within the volume of the electron wavefunction in a QD. We use a long rf pulse 
leading to selective (and complete) depolarization of the i-th isotope, thus enabling a selective 
measurement of the corresponding hyperfine shift (the change in the exciton spectral splitting) 
AEz,i- The AEz,% can be expressed via the nuclear polarization degree Pn,i of the i-th isotope as: 

AE z ,i = PiAiIP N:i , (S9) 

where pi is the relative concentration of that isotope and A{ is its hyperfine constant^. For 
InP, where the contribution of 115 In dominates, the hyperfine constant has been measured 
experimentally^: Ai\h Iri « 47 ^eV. Using this value, and neglecting variation of electron 
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density between gallium and indium sites, we can estimate the hyperfine constant for 69 Ga as 
^ 69 Ga = ^( 115 /n) ^lls^ f/l ~ ^1 P e ^i an d similarly for 71 Ga as A7i Ga 65 peV. We also take into 
account that the two gallium isotopes have natural abundances cre9 Ga ~ 0.6 and <J7i Ga ~ 0.4. By in- 
troducing the total gallium concentration p Ga we can write pe9 Ga = cr ^69 Ga ^p Ga , P^Ga = a ( n Ga)PGa- 
Since 115 In and 69 Ga have similar NMR frequencies, we can only measure their combined Over- 
hauser shift AE z &9 Ga + AE z iis In . The Overhauser shift for 71 Ga, AE z n Ga , is measured sep- 
arately. Finally, the nuclear spin polarization for each isotope Pjvi can be calculated using the 
Boltzman distribution (Eq. IS7|) and thus expressed in terms of polarization parameter (3 for a 
given spin /. Using Eq. [39] we can write the following system of equations for p Ga , pin and (5: 

9 3 

AE Z 69 Ga + AE Z 115 In = -p/ n A(115 /n )P A r(115 /n ) + -(T^9 Ga ^p Ga A^9 Ga )P N ^9 Gci) , 

3 

AE z n Ga = -0"(7i Ga )PGa-4(7i G . a )P A r j (7i Ga ), 

PGa + Pin = 1- (S10) 



For InP/GalnP QDs we measured the following values of the hyperfine shifts: AE z ii5j n + 
^E Z 69 Ga ps 120 fieV and AE z ji Ga » 8 ^eV. Solving Eq. ETO] we find p Ga rj 35%, p /n rj 65% 
implying significant penetration of the electron wavefunction into the GalnP barrier and/or dif- 
fusion of gallium into the dot. We also find (3 ~ 0.8 corresponding to the average electron spin 
of \(s z )\ ~ 0.2. For InGaAs quantum dots we assume the same values of hyperfine constants 
Ai. Using the measured shifts AE z usi n + AE Z 69 Ga ~ 56 peV and AE z n Ga ~ 18 peV, we find 
f3 ~ 0.8 and, as expected from peaks amplitudes of the NMR spectrum in Fig. [3^c), a much lower 
concentration of indium pj n ~ 20%. For both types of quantum dots we find very similar high 
degrees of optically pumped nuclear spin polarization: Pjv,( 116 /n) ~ 0.8 and -P/v,( 71 Ga) ~ 0-6- 

There are two possible scenarios resulting in relatively low average concentration of indium 
(pin ~ 20%) in InGaAs QDs revealed by the NMR measurements: (i) the nominal InAs QD 
contains significant amount of gallium due to diffusion during the growth process, or (ii) the 
dot itself consists mainly of indium but has a small size resulting in weak charge localization and 
significant penetration of electron wavefunction into the GaAs barrier. In both cases QDs will have 
a larger band-gap compared to the large-size InAs dots. This is indeed observed in the experiment: 
the wavelength of the ground state exciton recombination in the studied dots is A ~ 915 nm 
compared to A > 1000 nm in large indium- abundant dots. 
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S6. ESTIMATION OF STRAIN IN A QD 

Elastic strain is described by a second-rank symmetric tensor ey. The resulting electric field 
gradient (EFG) Vij can be related to eij via a fourth-rank tensor Sijki as s * 2 

Vij =^S ijk ie kU (i,j,k,l = x,y,z). (Sll) 
k,l 

In a crystal with cubic symmetry there are only 3 independent components denoted as Su, S12 
and £44. The symmetry relation S12 = —Su/2 is usually used to account for the zero trace of the 
Vij tensor, leaving only 2 independent components^ 2 -. Thus Eq. IS 1 II can be rewritten as 

Vu = Sn(eu - (e^ + e kk )/2), j ^ k 

Vij = 2Su^ij, i ^ j. (S12) 

It follows from Eq. IS12I that isotropic (hydrostatic) strain e/j = e xx + e yy + e zz produces no 
EFG, while biaxial (including uniaxial) and shear strains give rise to quadrupole shifts of NMR 
resonances. 

For uniaxial (e xx = e yy ) strain = e zz — (e xx + e yy )/2 along the direction of the applied magnetic 
field, the frequency shift of the 1/2 o 3/2 transition from the CT reads as vq = ^hupi-i) ' wnere 
Q is the quadrupole moment, h is the Planck constant. For the studied isotopes the quadrupole 
moments are: Q( 69 Ga) w 0.17, Q( 71 Ga) ra 0.10, Q( 75 As) ps 0.31 and Q( 115 In) ps 0.8 barn (1 barn 
= 10 -28 m 2 ). For 71 Ga in InGaAs dots the maximum vq ~2.5 MHz can be estimated from the 
width of the sidebands in the NMR spectrum in Fig. E^e). Using the value Su ~ 2.7 x 10 22 V/m 2 
measured for gallium in bulk GaAs^i 2 - we estimate |e&| ~ 6%. 

In InP dots we use NMR on 115 In satellite transitions to estimate strain magnitude. According to 
Eq. IS81 the amplitude of the "inverse" NMR signal from the satellite transition I' z o I' z +1 decreases 
with increasing spin \l' z \, leading to insufficient signal for large \I' Z \ and unreliable estimation of the 
maximum quadrupole shift. This is overcome in an additional experiment, where long broadband 
rf pulse (without the gap) centered at the frequency of the indium CT transition is used to erase 
polarization of indium. The total width of the excitation spectrum w exc is varied. We find that the 
magnitude of the erased nuclear polarization initially increases with w exc and saturates at a constant 
level at w exc ~ 20 MHz. This allows the maximum shift of the 7/2 9/2 transition frequency 
to be estimated as ~ 10 MHz, which is also equal to 4z^q. Using then Su ~ 5.9 x 10 22 V/m 2 for 
indium in InP we estimate the maximum strain as |e&| ~ 5%. 

Elastic strain also affects the CT frequency. Using vq ~2.5 MHz derived for 115 In we find that 
the shift of 115 In resonance according to Eq. [S4] can be as large as 1/2 ~ — 0-65 MHz at 
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vl ~ 49.7 MHz. For example, in the case of r/ = and angle ~ 54.7° corresponding to zero shift 
of ST bands the second-order shift of the CT is v_u 2 1/2 ~ — 0.25 MHz. On the other hand in 
experiment [Fig. E^b)] we observe shifts only on the order of ±50 kHz compared to the frequency 
vi n in unstrained InPS^. This suggests that the deviation between the strain axis and external 
field (characterized by 9) as well as non-axial symmetry of EFG (characterized by rj) are small. 
Thus, the most likely reason for the inhomogeneous distribution of ST shifts is the variation of 
vq within the dot volume caused by variation of the biaxial strain magnitude e^. In particular, 
non-zero amplitude of ST bands at CT frequency observed in Fig. [3^a) for 115 In can be explained 
by complete biaxial strain relaxation at the center of the do t ' S15 resulting in uq = 0. For 75 As 
nuclei large EFG can result not only from elastic strain but also from electric fields created by 
random substitution of gallium atoms by indium. This may explain the further broadening of 75 As 
CT transition. 

We note a significant difference in the asymmetry of the NMR spectra in Figs. [3](e,f): for 75 As 
the low (high) frequency ST band is enhanced for cr~( + ) pumping while for indium and gallium 
isotopes this asymmetry is reversed. The sign of the asymmetry is determined by the relative signs 
of Zeeman and quadrupole contributions to the nuclear Hamiltonian Eq. IS 1 L Since all isotopes 
of InGaAs have positive gyromagnetic ratios 7 and positive quadrupole moments Q, the opposite 
asymmetries of the spectra can be attributed to the opposite signs of the electric field gradients 
Vz'z' (see Eq. IS2j) experienced by nuclei of anions (As) and cations (Ga or In). 



S7. SPIN COHERENCE IN THE NUCLEAR SPIN ENSEMBLE SUBJECT TO 

INHOMOGENEOUS STRAIN 



This section details the experimental procedure and the model used for analysis of the de- 
pendence of the nuclear polarization dynamics in InP dots on the spacing between the modes in 
the broad-band rf excitation. The dynamics of the rf-induced nuclear spin polarization decay is 
measured using a pump-probe technique. Initially, nuclear spins are polarized optically. Then an 
rf-pulse of duration T r f is applied after which the nuclear spin polarization is probed optically by 
measuring the exciton Zeeman splitting E z . The spectrum of the rf excitation (without the gap) 
consists of a large number of spectral modes with equal spacing w m and with equal amplitudes B m 
changed as B m oc ^Jw m to keep the total power of the rf pulse constant (See Fig. HJ). The spectral 
range 44-47.5 MHz of the rf excitation is kept fixed and corresponds to satellite transitions of In 
nuclei [Fig. [31(a)] . The experimental dependences of the Overhauser shift on the rf excitation time 
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T r j are shown with symbols in Fig. IS5t a) for two different mode spacings w m of the rf-excitation 
spectrum. As seen, the increase of w m from 0.2 to 4 kHz results in a significant increase of the 
time required to equalize populations of the nuclear spin levels coupled by the rf field. We use ex- 
ponential fitting with a time constant 1/R r f to quantify the rate R r f of rf induced depolarization. 
The resulting dependence of R r f on w m is shown in Fig. H{b) an d is repeated with symbols in Fig. 

EKb). 

We now present in full the model that allows the study of coherence properties of nuclear spins 
using nuclear spin dynamic measurements. We start by noting that, in high external magnetic 
fields, nuclear spins are characterized by very long T\ times« exceeding 100 s, arising from sup- 
pressed energy relaxation leading to low probability of spin flips. Intrinsic linewidths of NMR 
transition correspond to much shorter times T<i determined by energy-conserving decoherence of 
transverse components of nuclear spins. This decoherence is induced by random dipole-dipole fields 
of other nuclei, and is characterized by the autocorrelation function which we choose in the form: 

1 a 

g(t) = - exp[-t/T 2 ] - exp -*/ (oT 2 , (S13) 

1 — a 1 — a 

where a < 1 is a dimensionless parameter. For a = this corresponds to an exponential corre- 
lation function with the relaxation time T 2 derived in the general theory of relaxation caused by 
fluctuations*^. The addition of the second term in Eq. IS 131 allows the condition </(0) = required 
by time reversal symmetry to be satisfied. The lineshape of the nuclear spin transition is found as 
the Fourier transform of g(t): 

~, \ ?Hl + a) /o U1 

9{V) l + (l + Q) 2 (2^) 2 r 2 2 + a 2 (27r^) 4 T 2 4 ' 1 ' 

As expected for a = this expression corresponds to a Lorentzian lineshape. 

For a nuclear spin transition between I' z and I' z + 1 states at a frequency i/ MC , a small amplitude 
(non-saturating) rf field will result in depolarization, which can be described by the differential 
equation for the population probabilities d(pii z+ i — pj/J/dt = —W x (pi' z +i — Pi')- For broad- 
band rf excitation consisting of discrete modes with frequency spacing w m , each inducing an equal 
magnetic field B m , depolarization rate W is given by: 

N 

W(v nuc ) = g(vnuc -vo-kx w m ), (S15) 

fc=0 

where the summation goes over all modes of the rf field with frequencies Uk = Vq + k X w m (vq is 
the frequency of the first spectral mode) . Since the total width of the rf band is much larger than 
w m or the nuclear spin transition width w nuc , the summation in Eq. IS15l can be extended to ±oo. 
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FIG. S5. (a) Experimentally measured nuclear spin polarization decay (symbols) for the rf-excitation mode 
spacings w m =0.2 kHz and 4 kHz. Nuclear polarization is calculated as a difference of exciton spectral 
splittings Ez measured after rf-pulse with duration T r f and before this pulse. Lines show model fitting, 
(b) Dependence of rf induced nuclear depolarization rate R r f on the mode spacing w m . Values extracted 
using exponential fitting of the experimental decay curves are shown with symbols. The dashed (solid) line 
shows results of calculations using the presented model for a Lorentzian (non-Lorentzian) lineshape a = 
(a « 0.18). 

For each individual pair of nuclear spin states, the rf field will induce an exponential decay of 
polarization. However, the QD contains a large number of nuclear spins with randomly distributed 
transition frequencies v nuc . Thus, in order to calculate the dynamics of the total nuclear polar- 
ization, we need to average over all possible values of v nuc . Since the spectrum of the rf band 
is a periodic function (i.e. the modes are equally spaced), such averaging can be done over one 
period. We also assume a uniform distribution of the nuclear transition frequencies i^nuc- 

Thus, 

the following expression is obtained for the time dependence of the difference between the current 
and initial nuclear spin polarization, describing the dynamics of rf-induced depolarization: 

[Wm / oo \ 

P N (t) - Pjv(0) = -Pjv(O) + Pjv(0) J o exp 1-tB g(u nuc -kx w m )\ du nuc , (S16) 

where B 2 oc B 2 m jw m is proportional to the spectral density of the total rf field power, which is 
kept constant in our experiment when w m is varied. 

Using Eq. IS16I we can fit nuclear spin decay curves measured for different values of w m . In 
Fig. 185( a) this is demonstrated for w m =Q.2 and 4 kHz. From the fitting of the whole set of decay 
curves measured with rf excitation in the spectral range 44-47.5 MHz for different w m we find 
good agreement with experiment for a nuclear spin dephasing time T2 ~ 2.5 ms and lineshape 
parameter a ~ 0.18. The same magnitudes of these parameters are obtained from fitting the data 
of a similar experiment but for an rf excitation band in resonance with a different section of the 
satellite transitions at 47.0-49.0 MHz. 
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We can now use the results of the fitting to analyze coherent dynamics of indium nuclear spins. 
We start by noting that for small w m all nuclear spin transitions are excited with nearly the same 
efficiency (similar to a white-noise excitation) resulting in an overall decay close to exponential. 
For larger spacing, the model given by Eq. IS 161 significantly deviates from the exponential decay: 
at the initial stage those transitions which are in resonance with the rf modes quickly become 
depolarized, whereas the decay is significantly slower for off-resonance transitions. For the analysis 
of the lineshape parameter a we characterize the rate of the decay given by Eq. IS 161 using the 
time ti/2 where nuclear polarization decreases to 1/2 of the initial value. The decay rate can be 
estimated as R r f = log(2)/t 1/ / 2 (in the limit of small w m it coincides with exponential decay rate). 
The dependence R r f(w m ) calculated in this way from Eq. IS 161 with parameters obtained from the 
fitting (T2 ~ 2.5 ms and a ~0.18) is shown in Fig. IS5f b) with a solid line. Calculation with the 
same Ti but for a=0 is shown with a dashed line. 

It can be seen that in the limit of small w m decay rate does not depend on a for both a=0 
and a=0.18. However, at larger mode spacing modeling using pure Lorentzian lineshape fails to 
describe the experiment (i.e. it is impossible to fit simultaneously all decay curves corresponding 
to different w m using Eq. IS 161 with a = 0). This deviation is described by the fourth order term 
oc z/ 4 in the denominator of Eq. IS 141 This term is responsible for much smaller depolarization rate 
at large detunings v 3> 1/(71" T2) in case of non-Lorentzian lineshape. At large mode spacing w m 
most nuclear spin transitions are excited only via their "wings" (i. e. at large v), and thus the 
effect of non-Lorentzian shape becomes more pronounced as seen in Fig. IS5f b). Strong deviation 
of the nuclear transition lineshape from Lorentzian function reveals non-exponential character of 
nuclear spin decoherence (described by Eq. IS 131 in our model) and demonstrates the potential of 
the presented technique for deeper understanding of the spin coherence of quadrupole nuclei. 
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